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Abstract. We proved uniqueness and instability of the symmetric subsonic- 
sonic flow solution of the compressible potential flow equation in a surface 
with convergent areas of cross— sections. Such a surface may be regarded as 
an approximation of a two— dimensional convergent nozzle in aerodynamics. 
Mathematically these are uniqueness and nonexistence results of a nonlinear 
degenerate elliptic equation with Bernoulli type boundary conditions. The 
proof depends on maximum principles and a generalized Hopf boundary point 
lemma which was proved in the paper. 



1. Introduction 

This paper is a continuation of our study on transonic flows in nozzles |f 71 118j . 
A fundamental difficulty to this problem is the lack of special solutions in a phys- 
ical nozzle. However, as motivated by |14j . it is rather easy to construct various 
physically interesting special solutions, such as those (supersonic-subsonic) tran- 
sonic shocks and (subsonic-supersonic) transonic flows observed in physical nozzles 
[4], for the steady compressible Euler equations in certain Riemannian manifolds 
|17[ 118] . Therefore these manifolds may be regarded as approximate nozzles and 
the study of the special flows in them would be of certain importance for the un- 
derstanding of flows in physical nozzles. In the previous works [Sj 1101 we and 
our collaborators have studied stability and uniqueness of transonic shocks for the 
potential flow equation or the complete Euler system, which involve free boundary 
problems of uniformly elliptic equations (for potential flow equation) or elliptic- 
hyperbolic composite system (for Euler system) . In [171 IE] we construct special 
transonic flow solution and show its stability for potential flow equation, which is a 
small perturbation result on nonlinear equations of mixed type. The present paper 
is devoted to the uniqueness and instability of subsonic-sonic flow in a convergent 
approximate nozzle in the class of all functions satisfying certain reasonable 
restrictions. These are global results on degenerate elliptic equations. We first for- 
mulate the problem. 

Let S^ be the standard unit circle in R^, and A4 the Riemannian manifold 
{{x,y) G (0, 1) X S^} with a metric G = dx ^ dx + n{x)'^dy dy. Here n{t) is a 
positive smooth function on [0, 1] satisfies: 



Date: September 14, 2009. 

2000 Mathematics Subject Classification. 35J70, 35B50, 76H05. 

Key words and phrases. Hopf lemma, degenerate elliptic operators, potential flow, surface, 
subsonic, sonic, nozzle. 



2 



PAN LIU AND HAIRONG YUAN 



{Hi) (i) n"{t) > 0; (ii) n'(<) < for t e (0, 1), n'{l) = 0. 

Without loss of generality, we also assume that n(l) — 1. We call A4 a convergent 
approximate nozzle. In the following, we take (x, y) £ M = (0, 1) x (0, 27r) as a 
coordinate system in M.. Note that on y = and y = 271 there should pose periodic 
conditions. 

The potential flow equation governing steady irrotational isentropic perfect gas 
flows in M is (c.f. [18]) 

+n{x)n'{x) (c^ + -^{d2(p)A di^ = 0. (1.1) 



n{x) 

(Here and hereafter we set di = dx, d2 = dy.) Direct computation yields that this 
equation is of elliptic type if the flow is subsonic (c^ > {diip)"^ + {82^)'^ / n{xY) , and 
degenerate if the flow is sonic (c^ = {dw)"^ + (92</?)^/n(a;)^). We may compute the 
speed of sound c by the Bernoulli's law 

c2 , {d,^f + {d2V?ln{xf Co ^^2) 



7-1 2 7-1 

with 7 > 1 and cq > two given constants. Let p, p be the pressure and density of 
the flow respectively. We use p — a,s the state function and then c = \f'Uf^- 

We call Sfe = {fc} X S^, fc = 0, 1, respectively the entry and exit of the approxi- 
mate nozzle JVi. Obviously dM. = Sq U Si. In the following we assume that 

{H2) 9i(y5 > on both Eq and Ei, 

which means exactly that the gas flows in M. on Eq and flows out M on Ei. 
We study subsonic-sonic flows in M. The boundary conditions are 

+ = 5]o, (1.3) 

{diV? + {d2^? = If^ on El, (1.4) 



where 60, &i are positive constants. Note that ()1.3p and ()1.4p are Bernoulli type 
conditions and the choice of hi shows that the flow is sonic at the exit. We choose 
60 to satisfy the algebraic equation 

br\co-^^hl)n{Qr-'^bl+' (1.5) 

and 60 < ''i so that the flow is subsonic at the entry. It is straightforward to show 
that there exists such a unique &o- We pose these nonlinear Bernoulli conditions 
rather than Dirichlet or Neumann conditions, since from physical point of view, 
prescribing density (pressure) at the entry and exit are more reasonable [3J H] . By 
Bernoulli's law (|1.2p . this consideration leads to (|1.3p and (|1.4p . 
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Now suppose the flow is symmetric, that is, depends only on x. Then ()1.1|) - ()1.4|) 
are reduced to this problem of u{x) ~ ^'{x): 

d^{npu) = 0, a;e[0,l], (1.6) 

^ + .e[0,l], (1.7) 

7 — 1 2 7 — 1 

u = ba, x = 0, (1.8) 

u^bi, x = l. (1.9) 

Under the above choice of bo, the velocity u can be easily solved and the obtained 
flow is subsonic in U Sq, sonic on Si, and always accelerates in A4 (i.e., u'{x) = 
(p"{x) > 0). Therefore we may choose a potential function (pi, = ftix) with ip'f^ — u 
being a special solution to problem (|l.ip - (|1.4p . 

Apart from stability, it arises also naturally the question that whether this special 
solution is unique in the large. That is, under suitable assumptions, whether any 
solution to (|l.ip - (|1.4p . modulo a constant, must be the (fb obtained above. The 
following results provide a positive answer to the uniqueness, but a negative one to 
the stability. 

Theorem 1.1. (i) The onlyC^{M) solution to problem (|l.ip - (|1.4p which satisfies 
{H2) and 

diiip > on Si (1.10) 

is the symmetric solution ipb modulo a constant. 

(ii) Under the same assumptions as in (i), but replacing (jl.3p by 

[div? + = B{y) on So, (1.11) 

with B{y) a positive function not identical to b^ and satisfying either B{y) > b^ or 
B{y) < bg, then there will be no any solution to the problem (jl.ip . (jl.2p . p. lip 
and HTM. 



We remark that there is no any extra requirements such as the flow to be sub- 
sonic in M to show the uniqueness. The only assumption ()1.10p means that the 
flow is accelerating at the exit. From the study of transonic flows [9], and the non- 
uniqueness results claimed by many authors based on numerical simulations (c.f. 
[15j). it seems that this assumption is necessary. Assertion (ii) shows particularly 
that the subsonic-sonic flow is not stable under perturbations of the density (or 
pressure, or speed) of the gas at the entry, provided the flow is accelerating at the 
exit, where it is sonic. 

The rest of this paper. Section 2, is devoted to the proof of Theorcni ll.il The 
key point is to prove a version of Hopf boundary point lemma applicable to points 
on characteristic degenerate boundaries for degenerate elliptic equations, which is 
Theorem 12.11 stated in Section 2.1. Although there are many impressive progresses 
in the study of degenerate elliptic equations and mixed type equations in gas dy- 
namics and other fields in these years (see [2 [U [51 [71 [HI [HI [HI [131 [IS] and references 
therein), this generalized Hopf lemma seems to be new. It captures the remarkable 
property that for subsonic-sonic flow or transonic flow, the potential flow equation 
behaves like the heat equation —diip-\- 822^ = near the sonic line (the line where 
the equation is degenerate), and hence the lower order term —diip is essential in 
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studying these degenerate elliptic or mixed type equations (c.f. [21 [18]). This ob- 
servation, obtained by studying special subsonic-sonic flows and transonic flows in 
an approximate nozzle, would be important for studying the flows in a physical 
nozzle. 



2. Proof of Main Theorems 

2.1. A generalized Hopf Lemma. We first present a version of Hopf lemma ap- 
plicable to degenerate elliptic operators and characteristic boundary points, which 
may be regarded as a generalization of results like Lemma 3.4 in 5 or Lemma 7.1.7 
in [16] and is of independent interests. 

Let D be a bounded domain of class in M" with boundary dD, and let 

n n 

L := ^ a*^'(x)9„- + ^ b'{x)d^ + c{x) (2.1) 

i^j — l i—1 

be a second order degenerate elliptic operator in D (that is, a'^^{x)^i^j > for 
all ? = ($!,■■■ i£.n) S M" and x € D). In what follows we assume that a^^x), 
6*(a;), c{x) are all bounded and continuous functions in D. Let P be a point on 
dD and assume, in a neighborhood U C K." of P, the boundary dD has the form 
{x = (xi, ■ ■ ■ ,Xn) G K" : Xn — f{xi, ■ ■ ■ ,Xn-i) — 0} with / a function, and 
U DD lies in {xn > f{xi, • ■ • , a;„)}. Then u — (— • ■ ■ , — g,f-^ , 1) is an interior 
normal vector to dD. We have the following result. 

Theorem 2.1. Let u be a C'^{D) n C^(i)) Junction and attains a strict minimum 
[maximum) at P £ dD. Also assume that, in a neighborhood U Cl D of P we have 
c = and 

n — l f. n—1 o2 r 

i—l ^J — 1 

Then if there holds 

Lu<0 (> 0) in D, (2.3) 

we have 

^ > (< 0) at P. (2.4) 

If c <Q, the same conclusion holds provided that u{P) < (> 0). 
Proof. First note that we may choose the following transformation 

2/1 = Zl - Xi{P), y2= X2- X2{P), • • ■ , Vn^ - f{xi, ■ ■ ■ ,X„_i) (2.5) 

to map P to O in y-coordinates and straighten dD locally to y„ = near the origin 
O. Denote the image of U to be V. 
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It follows that the operator L can be written in the form 

dxi 




k,l = l \i,i = l 

where a"(,) := E^., ""feS; ""t" ^= E?- '"fe + Elj.i ""sgi". 

Hence by the assumption (|2.2[) we have 




n— 1 ^ n — 1 

> 0. (2.7) 
Now we introduce an auxiliary function h of the form 

n 

h{y) = -Y^yl + liyn, (2.8) 
1=1 

where is a positive constant yet to be determined. A direct calculation gives 

n 

M(0) = -2^a"(0) +^/3"(0) (2.9) 

i=l 

and 

|^(O)=/^>0. (2.10) 

Note that we have /3"(0) > 0. Hence fi may be chosen large enough so that 
Lh{0) > 0. 

Let us pick up a small rectangle Di := {—di,di) x (— d2,rf2) x ••■ x (0, c) C V, 
with O G and Lh > through out the domain Di. Corresponding to the 
function h there is an open set B — {y € R" : h{y) > 0}. One sees from the 
construction of h that K :— Di O B is non-empty, and the boundary of K can 
be divided into two parts: Ki := dK n dDi, K2 := dK n dB, with O £ K2. 
Now since u — u(0) > on /Ci, there is a constant e > for which the function 
w :— {u — u{0)) — eh > Q on Ki. This inequality is also satisfied on K2 where 
h = 0. Thus we have Lw — Lu — cu{0) — eLh < —eLh < in K, and w > on dK. 
The weak minimum principle (see the first paragraph of the proof of Theorem 3.1 
in [5]; it is here we need the matrix (a*^ ), hence (q:'^')j to be semi-positive definite) 
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now implies that w > in K. Taking the interior normal derivative at the point 
O, where w{0) — 0, we have -§^{0) > 0, which impHes 

|i,0,>.^,0)...>„^ ,2.1, 

Therefore we get, in the original x-coordinates, 



> (2.12) 

du 



as required, since there should hold §^{0) = for i = 1, • • ■ , n — 1 due to the fact 



that u attains a minimum at O. 

The case that u attains a strict maximum at P G dD can be proved by applying 
the above proved result to —u. For c < 0, the proof is similar. □ 

Remark 2.1. As pointed out in [5] (p. 34), (|2.4|) can be replaced by 

^.^ K.)-.(P)| ^^ (2.13) 
x^p \x — F\ 

where the angle between the vector x — P and the inner normal vector of dD at P 
is less than 7r/2 — 5 for some fixed positive number 5. We require u S C^(Z)) is in 
essence just to obtain that the tangential derivatives of u along dD are zero at P, 
which is used to derive ()2.12|) . 

2.2. Proof of Theorem 11.11 We first prove part (i) of Theorem ll.il 

Suppose ip \s a. C'^ solution to problem (|l.ip - (|1.4p . Then by (jl.ip . it is straight- 
forward to check that tjj ^ iph — (p satisfies the following equation 

2 

J2 a'H.^b)d^ji> + b^dii> + ^"82^ = 0, (2.14) 



where 



a^\^b) = n{xf{cl - [di^bf), (2.15) 

a^^iifib) = a^^i^b) = -diVbd2Vb = 0, (2.16) 

a^^^b) ^ cl - -^{d2Vbf ^ cl (2.17) 
nyx)'^ 



with Cb — y cq — ^^^(3i(^h)^ the sonic speed corresponding to ifb, and 
- - (^^^n{x)^dn^ + ^a22(^) di{^ + ^b) 

+n{x)n'{x) (^cl + ^^^di^di{^ + V'fc)) , (2-18) 

,7—1 7+11 

h = — -— 92v9ii<^ + 2ditpdi2(p H rT7^2'^922<p 

2 2 n{xY 

+ ^^d,^d2^. (2.19) 
2 n(x) 

Recalling that ipb is a subsonic-sonic flow, so equation (j2.14p is a linear degenerate 
elliptic equation of ip: It is strictly elliptic in 7W U So and degenerate on Si. 
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The boundary conditions of 'ip are 

lo-Dtp^Q on So, (2.20) 
h-Dtp^Q on El, (2.21) 

where D-ip — (diip, 82^^) and 

By {H2) and the fact diipi, > 0, p.20p and ()2.2ip are both obhque derivative con- 
chtions. 

Now we apply maximum principles to the problem (|2.14p . (|2.20l) and (|2.2ip . If 
ip is not a constant, by strong maximum principle which is valid in any compact 
subset of A4, the maximum and minimum of tjj can only be achieved on Eq or 
El. By the classical Hopf boundary point lemma (Lemma 3.4 in [5J), it is only 
possible to achieve the extremum on Ei, i.e., the boundary where the equation is 
degenerate. 

To show Ip is a constant, by Theorem 1 2 . 1 1 and an argument of contradiction, we 
just need to show the validity of ()2.2I) . In our situation, it is exactly 

6^ < on El. (2.22) 

The calculation is straightforward and we omit it (c.f. (|2.7I) and also note that 
y„ = -{x - 1) now). 

Since diipb > 0, diip > and n'(l) = 0,n(l) = 1 by {Hi){H2), ^(TM is guaran- 
teed by 

2±}_d^^^+2^d22^>Q on El. (2.23) 

We show this is true under the assumption (|1.10p . 
In fact, by the equation Ijl.ip . we obtain that on Ei, 

c^l\ip = {diipfdiiLp + 2diLpd2^di2ip + (d2^fd22^- (2.24) 

Differentiating the boundary condition (|1.4p with respect to y and multiplying 82^^ 
it follows 

diipd2ipdi2ip + {d2vfd22V = 0. (2.25) 

Substituting this in ()2.24p and note that = [diip)"^ + {d2^Y , then 

(c^ + {d2ipf)l^^ = c^^ii^ > 0. (2.26) 

So Av3 > on El and (P?^ follows. (Note that c = 61 > on Ei.) This finishes 
the proof of part (i) of Theorem 11.11 

We continue to prove part (ii) of Theorem ll.il The only difference is the bound- 
ary condition (|2.20p . Now it should be written as 

l^-DiP^hl-B{y) on Eq. (2.27) 

We consider the case B > bg. Suppose tp is not a constant. Then by the classical 
Hopf lemma, the minimum of ip can only be achieved on Ei. But this is contra- 
dictory to (|2.22p and Theorem 12.11 So ^p must be a constant. Hence on Eg there 
should hold 69 = B{y). However, we have assumed that B is not identical to fcg, 
so the only possibility is that there is no any solution of the problem (jl.ip . 
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(|1.2|) , and ()1.4p . The case B < can be proved similarly by considering the 

maximum of "0. 

This finished the proof of Theorem 11.11 

Remark 2.2. We see the equation (|2.14p is reduced to 

cld22'4' + h^dii! + h^d2il^ = (2.28) 

on El. So (j2.22p implies that this equation is similar to a heat equation. This type 
of degeneracy also occurs, for example, in shock reflection phenomena (c.f. 
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